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Statistical mechanics of ID multivalent Coulomb gas may be mapped onto non-Hermitian quan- 
tum mechanics. We use this example to develop instanton calculus on Riemann surfaces. Borrowing 
from the formalism developed in the context of Seiberg-Witten duality, we treat momentum and 
coordinate as complex variables. Constant energy manifolds are given by Riemann surfaces of genus 
g > 1. The actions along principal cycles on these surfaces obey ODE in the moduli space of the 
Riemann surface known as Picard-Fuchs equation. We derive and solve Picard-Fuchs equations for 
Coulomb gases of various charge content. Analysis of monodromies of these solutions around their 
singular points yields semiclassical spectra as well as instanton effects such as Bloch bandwidth. 
Both are shown to be in perfect agreement with numerical simulations. 



I. INTRODUCTION 

One of the very last works of Anatoliy Larkin^ was 
devoted to transport through ion channels of biologi- 
cal membranes. An ion channel may be roughly viewed 
as a cylindrical water-filled tube surrounded by a lipid 
membrane. Its typical radius a ~ 6 A is much smaller 
than its length L w 120 A The important observa- 
tion with far reaching consequences, made in Ref. pQ, 
is that the dielectric constant of water e wate r ~ 80 is 
significantly larger than that of the surrounding lipid 
membrane eii p id ~ 2. This defines a new length scale 

f ~ Q>yJ ^water/ ^lipid m(e wate r/ eiipid) ~ 14 0A over which 
the electric field stays inside the channel and does not 
escape into the surrounding media. Since £ > L, the 
ions inside the channel interact essentially through the 
ID Coulomb potential U(x\ — x^) ~ eE§\x\ — x<i\, where 
Eq = 2e/a 2 e wate r is a discontinuity of the electric field 
created by a unit charge. This fact dictates a significant 
energy barrier U(L/4) w 4/c#T room for moving a single 
ion through the channel. If indeed present, such a bar- 
rier would essentially impede ion transport, preventing 
the channel from performing its biological functions. 

Nature removes such Coulomb blocking by screening. 
A moving ion is screened either by mobile ions of disso- 
ciated saltP, or by immobilized charged radicals attached 
to the walls of the channePJE Nevertheless, due to the 
peculiar nature of the long-range ID Coulomb potential, 
the transport barrier proportional to the channel length 
L is always present. Its magnitude, though, is typically 
suppressed 1 down to about &BT room , allowing for a rela- 
tively unimpeded transport of ions. These considerations 
call for development of a transport theory of ID Coulomb 
gases. Following the celebrated mapping of ID statisti- 
cal mechanics onto an effective quantum mechanics, pio- 
neered by Edwards and LenarcP^ and Vaks, Larkin and 
Pikin 11 , reference [1 mapped the problem onto quantum 
mechanics of a cosine potential (we briefly review this 
mapping in Sec. [il]). The ground state energy of such 



quantum mechanics is exactly the equilibrium pressure in 
the Coulomb plasma. Moreover the width of the lowest 
Bloch band is a specific energy barrier for ion transport 
through the channel. 

It is instructive to notice that 2a cos = tx(e %d + e~ z0 ) 
potential describes a mixture of positive, e* 6 *, and nega- 
tive, e~ %e ', monovalent ions with concentration a. One 
may also consider a situation when the channel is filled 
with a solution of dissociated multivalent salt, such as e.g. 
divalent CaCl2 or trivalent AICI3. In these cases the cor- 
responding ID statistical mechanics is mapped onto the 
quantum problem with a non-Hermit ian p otential such 
as a(\e 2ie + e~ ie ) or a(\e 3i0 + e~ ie ) 2 10 . The present 
paper is devoted to efficient mathematical methods of 
treating non-Hermitian quantum mechanics of this sort. 

Our particular focus here is on a semiclassical treat- 
ment, applicable in the regime of sufficiently large salt 
concentration a. In its framework the energy spec- 
trum (thus the pressure) is determined by the Bohr- 
Sommerfeld quantization condition for the action of clas- 
sical periodic orbits. On the other hand, the bandwidth 
(and thus the transport barrier) is given by the expo- 
nentiated action accumulated on the instanton trajec- 
tory, running through the classically forbidden part of 
the phase space. The traditional techniques of Hermi- 
tian quantum mechanics call for finding classical and in- 
stanton trajectories by solving equations of motion in 
real and imaginary time and evaluating corresponding 
actions. This route can't be straightforwardly applied 
to non-Hermitian quantum problems arising in the con- 
text of multi-valent Coulomb gases. Even leaving aside 
the technical difficulties of solving complex equations of 
motion, there are conceptual difficulties with identifying 
periodic orbits as well as the meaning of classically al- 
lowed vs. forbidden regions and with the imaginary time 
procedure. 

In this paper we borrow from the algebraic topology 
methods developed in the past decades in the context 
of the Seiberg-Witten solutiorf^^ and its applications 
to integrable systems^E^ (and many follow-up contri- 
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butions). The central idea is to consider both coordi- 
nate and corresponding canonical momentum p as com- 
plex variables. This leads to /b?/r- dimensional (4D) phase 
space. Then (complex) energy conservation restricts the 
trajectories to live on 2D Riemann surfaces embedded 
into 4D phase space. The dynamics of the system are es- 
sentially determined by the topology, i.e. genus g, of such 
Riemann surfaces. We show that e.g. mono- and divalent 
gases are described by tori, while trivalent and 4-valent 
lead to genus-2 surfaces, etc. The Cauchy theorem and 
the resulting freedom to deform the integration contour 
in the complex space allows to avoid finding specific so- 
lutions of the equations of motion. Instead one identifies 
the homology cycles on the Riemann surface and finds 
the corresponding action integrals, which depend only on 
the topology of the cycles and not on their specific shape. 
For example, the cosine potential of monovalent gas leads 
to a torus, which obviously has two topologically distinct 
cycles, Fig. [6] The two turn out to be related to clas- 
sical and instanton actions correspondingly. The genus 
g > 1 Riemann surfaces admit 2g topologically distinct 
cycles. Below we identify and explain the meaning of the 
corresponding action integrals. 

The shape of the specific Riemann surface depends on 
the parameters of the problem, e.g. salt concentration 
a in our case. Such parameters are called moduli of the 
Riemann surface. It turns out that the action integrals, 
being functions of the moduli, satisfy closed ordinary dif- 
ferential equation (ODE) of the order 2g, known as the 
Picard-Fuchs equation. The actions may be found as 
solutions of this ODE in the moduli space, rather than 
performing integrations over cycles on the surface. Be- 
low we derive and solve Picard-Fuchs equations for sev- 
eral (positive, negative) ionic charge combinations, such 
as genus g = 1 cases (1,1), (2, 1) and genus g = 2 cases 
(3, 1), (3, 2), (4, 1). We then discuss how to connect the 
principal classical actions with the spectra of the corre- 
sponding quantum problem. The key observation is that 
in the moduli space the actions exhibit a few isolated 
branching points. Going around such a branching point 
transforms the actions into their linear combinations - 
the Sp(2g, Z) monodromy transformation. The invari- 
ance of quantum observables under monodromy transfor- 
mations dictates Bohr-Sommerfeld quantization for one 
of the principal classical actions. The remaining actions 
may be identified with the instanton processes, related 
to e.g. Bloch bandwidth. 

Statistical mechanics of ID Coulomb gases may seem 
to be an isolated problem, not worthy of developing an 
extensive mathematical apparatus. Our goal here is to 
use it as a test-drive example, grounded into a well- 
posed physics problem, to develop a machinery applica- 
ble in other setups. Recently the so-called VT symmet- 
ric non-Hermitian quantum mechanics attracted a lot of 
attention for its application in active optics^ anc [ open 
quantum systems^, as well as in the description of an- 
tiferromagnetic lattices^ and calculating energy states 
in larger molecules 22 . Our examples also belong to the 



class of VT symmetric problems. It seems likely that the 
methods developed here may be applied to advance ana- 
lytical understanding of a broader class of VT symmet- 
ric quantum mechanics. Another context, where com- 
plexified quantum mechanics was proven to be ex tremely 
useful, is dynamics of large molecular spin J 23|24 4 Indeed 
functional integral representation of the spin dynamics 
leads naturally to the Hamiltonian formulation, where 
the projective coordinates (z, z) on the sphere play the 
role of the canonical pair 2 -^. It was realizecP^ 2 ^ that to 
find instanton trajectories one has to consider z and z 
as independent complex variables, thus expanding the 
dynamics into 4D phase space. The Riemann geometry 
methods seem to be well-suited to advance this subject 
as well. 

This paper is organized as follows: in section[ll]we out- 
line the relation between ID multivalent Coulomb gases 
and non-Hermitian quantum mechanics and discuss gen- 
eral symmetries of the latter. In section III we summarize 



major numerical observations regarding complex spectra 
and band-structure for the family of Hamiltonians con- 
sidered here. In section |IV| we illustrate the machinery 
of algebraic geometry on Riemann surfaces for the famil- 
iar Hermit ian cosine potential quantum mechanics, which 
corresponds to the monovalent (1,1) gas. Here we intro- 
duce complexified phase space and Riemann torii of con- 
stant energy; we then derive, solve and analyze solutions 
of the Picard-Fuchs equations. In section |V| we apply the 
developed methods for the divalent (2, 1) Coulomb gas, 
which is also described by a genus- 1 torus. In section fVT 
we extend the method for genus-2 example of trivalent 
(3, 1) gas, which exhibits some qualitatively new features. 
The (3, 2) and (4, 1) gases are briefly tackled in section 
|VII| In section |VIII| we outline connections to Seiberg- 
Witten theory. We conclude with a brief discussion of 
the results in section llXl 



II. MAPPING OF COULOMB GASES ONTO 
QUANTUM MECHANICS 

Consider a ID gas of cations with charge n\e and an- 
ions with charge — ri2e, where (711,712) are positive in- 
tegers. By Gauss's theorem, the electric field at a dis- 
tance x larger than the radius of the channel a from 
a unit charge is Eq = 2e/a 2 e water . At the location of 
a charge ni^ the electric field exhibits a discontinuity 
±2^o n i,2- Since all charges are integers the field is con- 
served modulo 2Eo along the channel. This allows us 
to define the order parameteiE^ q = E(x)( mod 2Eq), 
which acts like an effective boundary charge ±q at the 
two ends of the channel. The Poisson equation in ID 
reads V 2 </> = — 2Eq5{x) 1 leading to ID Coulomb poten- 
tial (j)(x) = — E \x\. The potential energy of the gas is 
thus 



U = 



eE 
2 



1,3 



(i) 
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where <jj is the charge rt\ or — n 2 of an ion at the position 
Xj and we omit the ±q boundary charges for brevity. Our 
goal is to evaluate the grand canonical partition function 
of the gas in the channel of length L 



Zr = 



Nx,N 2 =0 1 Z i=l J ° j = l 



-U/k B T 



(2) 

where fi^ are fugacities of the two charge species. 
One can now introduce the charge density employ- 
ing a delta-function S[p(x) — ^2j&jS(x — Xj)). The 
delta-function is elevated in the exponent with the 
help of the auxiliary field 0{x). This procedure de- 
couples all Xj integrals- 1 -, bringing them to the form 
J2 N [f fdxe ia °W] N /N\ = exp{f fdxe ia6 W}. The in- 
teraction potential ([I]), being inverse of the ID Laplace 
operator, leads to exp{(T/e£?o) JdxQd^O}. As a result 
the partition function Q is identically written as the 
Feynman path integral, in an "imaginary time" x, for 
the quantum mechanics with the Hamiltonian 



H = {ids - q) 2 - (a ie ini6 + a 2 e- in29 ) , 



(3) 



where ai^ — f 1,2k bT / eE<) are dimensionless ion concen- 
trations. Such Feynman integral is the expectation value 
of the evolution operator during "time" L, leading to 



Z L 



Xe 



■ Jq dx H 



\(q\m)\* 



»(<z) 



(4) 

where X stands for x-ordered exponent. Here e m (q) are 
eigenvalues of the effective Hamiltonian H and \m) = 
ipm(6) are its eigenvectors in the Hilbert space of periodic 
functions ip m (0) = VVn(# + 27r), and finally the matrix 
elements are (q\m) = j^dOe- iqe ^ m (e). The boundary 
charge q plays the role of the Bloch quasi- momentum 
and the spectrum is obviously periodic in q with the unit 
period (reflecting the fact that the integer part of the 
boundary charge may be screened by mobile ions and 
thus inconsequential). 

The pressure of the Coulomb gas is its free energy per 
unit length 



dlnZ L l- 
dL 



-eE e (q) , 



(5) 



where eo(q) is the eigenvalue with the smallest real part. 
In equilibrium the system minimizes its free energy by 
choosing an appropriate boundary charge q. In all cases 
considered below the minimum appears to be a non- 
polarized state of the channel, i.e. q = (see how- 
ever Refs. [2 for exceptions to this rule). Adiabatic 
charge transfer through the channel is associated with 
the boundary charge q sweeping through its full period. 
As a result, the (free) energy barrier for ion transport is 



U = eE LA , 



(6) 



where Aq is the width of the lowest Bloch band. There- 
fore the ground state energy and the width of the lowest 



Bloch band of the Hamiltonian (|3| determine thermody- 
namic and transport properties of the (711,712) Coulomb 
gas. The rest of this paper is devoted to a semiclassical 
theory of the spectral properties of such Hamiltonians. 
We start by discussing some general symmetries of the 
non-Hermitian Hamiltonian (|3|. 



A. VT Symmetry 

Although the Hamiltonian Q is non-Hermitian for 
^1 7^ ^2, it obeys VT- symmetry 26 27 . Here the parity 
operator V acts as 6 — ^ — while the time-reversal op- 
erator T works as complex conjugation i — » —i. Clearly 
the two operations combined leave the Hamiltonian (|3| 
unchanged. One may prove^ZL^ tn a t all eigenvalues of 
T^T-symmetric Hamiltonians are either real, or appear 
in complex conjugated pairs. As shown below for posi- 
tive values of concentrations > the lowest energy 
band eo(q) is entirely real, ensuring the positivity of the 
partition function. The higher bands e m (q) are in gen- 
eral complex. It is interesting to note that, for unphysical 
negative concentrations ai^ < 0, already the lowest band 
eo(g) is complex, making the free energy ill-defined. 



B. Isospectrality 

The spectrum of the Hamiltonian (|3| is invariant under 
shift of the coordinate — » 0+9o, where #0 is an arbitrary 
complex number. Upon such transformation (preserv- 
ing the periodic boundary conditions) the dimensionless 
concentrations ai 5 2 renormalize as ol\ —> aie %nie ° and 
0L2 — >• a2e -m2<9 °. Notice that the combination a^a^ 1 re- 
mains invariant. From here one concludes that the family 
of Hamiltonians ([3| with 

a™ 2 a™ 1 = const (7) 

are is spectral 10 . Thus without loss of generality, one 
may pick one representative from each isospectral fam- 
ily. It is convenient to choose such a representative to 
manifestly enforce charge neutrality in the bulk reser- 
voirs. To this end one takes a^ni = a2U2 = a, which 
brings the Hamiltonian (|3| to the form 



H = a 



P 2 ~ I — e' 
ni 



n 2 



(8) 



where we have defined the momentum operator as 

p = a- 1 / 2 (-id e + q); [0,p] = ia^ 2 . (9) 

The commutation relation shows that a~ x l 2 plays the 
role of the effective Planck constant. With the help of 
the isospectrality condition ([7]), one may always choose 
a proper a such that the spectrum of Hamiltonian ([8| is 
identical with that of a Hamiltonian with arbitrary a\^- 
The physical reason for this symmetry is that the inte- 
rior region of the long channel always preserves charge 
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neutrality, allowing the edge regions to screen charge im- 
balance of the reservoirs. Therefore, irrespective of the 
relative fugacities of cations and anions in the reservoirs, 
the thermodynamics of the long channel are equivalent 
to the one in contact with neutral reservoirs with an 
appropriate salt concentration a. Hereafter we restrict 
ourselves to the neutral Hamiltonian (|8| with the single 
parameter <x 



III. NUMERICAL ANALYSIS 

In this section we discuss numerical simulation of the 
spectrum of the Hamiltonian ([sj) - We focus on unequal 
charges n\ ^ 77,2, since the case of rt\ = re duces 
to the well-known Hermitian cosine potential 29 30 . For 
unequal charges the Hamiltonian is non-Hermitian but 
PT- symmetric, allowing for complex eigenvalues which 
appear in conjugated 

Since the Hamiltonian H acts in the Hilbert space of 
periodic functions, one may choose the complete basis 
in the form {e* m6> } mG z- In this basis the Hamiltonian is 
represented by an infinite size real matrhP 



Hrr 



■■(m-q) 2 6 n 



Tlx 



n 2 



(10) 

The boundary charge q plays the role of quasi-momentum 
residing in the Brillouin zone q G [— ^, \]. To numer- 
ically calculate the energy spectrum e m (q) we truncate 
the matrix at a large cutoff, after checking that a further 
increase in the matrix size does not change the low-energy 
spectrum. We left the boundary conditions "open", i.e. 
did not change the matrix elements near the cutoff, after 
verifying that different boundary conditions don't affect 
the result. It is easy to see that the matrix size should be 
^> yfot to accurately represent the low-energy spectrum. 
As an illustration we show the Hamiltonian cut to a 5 x 5 
matrix for divalent (2, 1) gas: 



/(-2-9) 2 
—a 



V 



-f 

(-l-q) 2 -f 

-a (0-q) 2 -f 

-a (l-q) 2 

-a (2-q) 2 ) 



For reasons which will become apparent below, it is 
convenient to present the spectrum e on the complex 
plane of the normalized energy u defined as 



n x n 2 e 
Ti\ + ri2 a 



(ii) 



For the divalent (2, 1) gas u = 2e/3a and the correspond- 
ing spectra are shown in Fig. [I] The spectrum consists 
of a sequence of complex Bloch bands. The number of 
narrow bands within the unit circle \u\ = 1 scales as yfa. 
They form three branches which terminate at u = —1 and 
u = e ±l7T / 3 and approximately line up along the lines con- 
necting the termination points with the point u = 1. We 



(a) a — 0.5 



(b) a = 1 



(c) a = 2 



(d) a = 200 

FIG. 1: (Color online) Complex plane of normalized energy 
u = 2e m (q)/Sa for (2, 1) gas. The color corresponds to 
different values of quasimomentum q; blue stands for q = 
and red for q = ±1/2. The dotted circle is \u\ = 1, the 
dashed lines connect between u — 1 and u — e ±l7r ^ 3 , 
indicating positions of the narrow complex bands in the 
limit of large a. 



shall discuss the corresponding bandwidths below. Out- 
side the unit circle the bands are wide and centered near 
the positive real axis of energy. 

Figure [2] shows the band structure in the first Brillouin 
zone \q\ < 1/2 for a = 1. Notice that the lowest Bloch 
band is purely real (this is always the case for a > 0), 
ensuring positive partition function Q and real pressure 
([5]). The next two bands are complex. For \q\ < q c « 0.36 
they exhibit opposite imaginary parts (not shown), but 
turn real at \q\ > q c . The next two bands are real, cf. 
Fig. [TJd. The higher bands form an alternating sequence 
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FIG. 2: (Color online) Band structure for (2, 1) gas with 
a — 1, cf. Fig. [T]d, vs. boundary charge (quasi- 
momentum) q. For the complex bands the real part of 
€m(q) is shown in dashed blue. 



(a) (2, 1) 



(b) (3, 1) 



of two real and two complex bands. For larger values of 
a there is a sequence of entirely complex narrow bands, 
cf. Fig.fltL 

FigurejS] shows normalized spectra for several different 
combinations of charges on the complex energy plane of 
u, Eq. (11), at large concentration a = 200. One may 
notice odd number ri\ + ri2 or m + ri2 — 1 of spectral 
sequences, consisting of order y/a exponentially narrow 
bands, seen as points. The central sequence goes along 
the real axis terminating at the bottom of the spectrum 
near u = —1. The other appear in conjugated pairs 
terminating near the roots of unity u = — (l) 1 /( ni + n2 ). 
Close to the termination points the band sequences align 
along the lines pointing towards u = 1. Further away 
from the termination points they deviate from these lines 
and may coalesce. 

Although thermodynamics and transport properties of 
the Coulomb gases are merely determined by the lowest 
band eo(q), below we address the wider spectral prop- 
erties of Hamiltonians ([8]), presented in Figs. [T] — [3j To 
this end we develop a semiclassical theory which is best 
suited for the description of exponentially narrow bands 
present at large concentration a > 1. 



IV. MONOVALENT (1,1) GAS 

To introduce the methods, we first develop a semiclas- 
sical spectral theory for the Hermitian Hamiltonian ([8|, 
([9| with ni = ri2 — 1. To this end we look for wavefunc- 
tions in the form ip — e l ° 1/2s , where S is an action for 
the classical problem with the normalized Hamiltonian 



(c) (4, 1) 



(d) (3,2) 

FIG. 3: (Color online) Complex plane of normalized energy 
u, Eq. |TT| ), for a = 200 and various valences (711,712). The 
dotted circle is \u\ = 1, the dashed lines connect spectrum 
termination points u — — (l) 1 /( ni + n2 ) anc [ u = 1, 
indicating positions of narrow complex bands. 



where u = e/(2a), so u = =pl correspond to the bottom 
(top) of the cosine potential. The semiclassical calcu- 
lations require knowledge of the action integrals. Our 
approach to such integrals is based on complex algebraic 
geometry. First, let z = e %d and consider (z,p) as com- 
plex variables. Since p(z) resides on the constant energy 
hypersurface 



2u = p 2 



(13) 



2u = p 2 - 2 cos , 



(12) 



we have a family of complex algebraic curves 

S u : F{p, z) = p 2 z - (z 2 + 2uz + 1) = (14) 
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parameterized by u. For u ^ =f1 it can be checked that 
(dT/dz, dT /dp) does not vanish on £ n , so each S u is non- 
singular. Then J-*(p, z) implicitly defines a locally holo- 
morphic map p = p(z). The exceptions to this occur at 
z = 0, oo, z±, where 



z± 



(15) 



are the roots of p 2 = (i.e. classical turning points). In 
a vicinity of these four branching points p(z) behaves as 



p ~ z 



-1/2 



p~z x '\ 
p~(z- z±) 1/2 , 



(z ~ 0) 
(z ~ oo) 
(z ~ z±) 



(16) 
(17) 
(18) 



respectively, i.e. p(z) is locally double- valued. (Note that 
we have added a point at z = oo to the complex plane, 
thereby rendering it compact and topologically equiva- 
lent to a Riemann sphere, Fig.|4|. To make sense of this 
double- valuedness, we first introduce two cuts between 
the four branching points. For convenience we have cho- 
sen to do so between 0, oo and the turning points z±. 
Upon this cut domain, p(z) is locally holomorphic. 




OC 




(a) 



FIG. 4: (a) Complex z-plane with two cuts, (b) It 
compact ifies to Riemann sphere with two cuts. 



FIG. 5: Construction of Riemann surface of genus 1. Two 
Riemann spheres with two cuts each are deformed into 
tubes to make the gluing in the final step more clear. 




We then introduce a second sheet of the z-plane and 
the corresponding Riemann sphere, cut in the same way 
as the first. We then analytically continue p(z) on the 
first sheet across the cuts onto the second sheet. If p(z) 
is analytically continued across the branch cut again, we 
arrive back on the first sphere where we started. In this 
way, we obtain p(z) as a locally holomorphic function, 
whose domain is a doubly-branched cover of the Riemann 
sphere. Furthermore, suppose we open up the branch 
cuts, keeping track where on the other branch p(z) will 
be, if we cross one side of a cut. Identifying these edges 
one obtains a torus as in Fig. [5] (where the arrows are 
used to signify the glued together edges) . Thus the com- 
plex algebraic curve S u can be understood as a compact 
Riemann surface of genus g = 1 (generically, every com- 
pact Riemann surface is topologically a sphere with some 
number of handles g, called the genus of the surface). 

In the exceptional points u = =pl the two turning 
points collide (z+ = Z- = ±1) and the branch cut be- 
tween them collapses. The Riemann surface degenerates 
into a sphere with two points identified, a singular sur- 
face of genus 0. This coincides with one of the loops of 
the torus becoming contractible to a point, Fig. [6] 



FIG. 6: Riemann surface of genus- 1 with two basic cycles So 
and Si on it. In the limit u —> =pl the torus degenerates 
into a singular surface. This coincides with the loop So 
(but not Si) becoming contractible to a point. 



A. Integration and topology on torus 

The action integrals can be understood as S = §^ A 
over classical trajectories, where 

\{u) = p{9) d0 = p(z) - = {Z2 + 2UZ + 1)V2 dz (19) 

IZ IZ^' Z 

is the action 1-form which is holomorphic on the torus. 
To visualize the relevant trajectories we momentarily re- 
turn to and consider it as complex. In this representa- 
tion one has square-root branch cuts along the real axis, 
connecting the classical turning points. The action in- 
tegrals run just above or below the real axis in between 
the turning points. Combining them into closed cycles, 
one can push these cycles off the real axis and away from 
the turning points without altering the action integrals 
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(by Cauchy theorem). The two deformed cycles, shown 
in Fig. [7| are hereafter called 70 and 71 . 




FIG. 7: The classically allowed (forbidden) region at energy 
2u are shown by the solid (dashed) gray line. A classical 
(instanton) periodic orbit, in the complex #-plane, leads to 
70(71) cycles. 



Translating these two cycles to the complex z-plane 
yields the contours of Fig. [8j Notice that these are in- 
deed cycles (i.e. closed contours) owing to the crossing of 
branch cuts. On the Riemann surface both wind around 
the torus. For this reason, the integrals Sj(u) = § A are 

known as periods of £ u with respect to \{u). Cauchy's 
theorem is also valid on the Riemann surface since the 



action form (19) is holomorphic there by construction of 



the torus. This has a major implication. Consider the 
cycles So, Si as defined in Fig. [6j Any closed cycle on 
the torus (after appropriate deformation) can be decom- 
posed into a superposition of an integer number of these 



two basic cycles. For example, the cycles 70 and 71 are 
7o=^o, 71 = 25! -S . (20) 

This is evident if one examines the manner in which these 
cycles encircle around the torus. Formally, the basic cy- 
cles generate the first homology group of the torus (since 
cycles which are alike in this manner are homologous). 

One can also consider the first co homology group of the 
torus, comprised of all holomorphic 1-forms on the Rie- 
mann surface modulo exact 1-forms (the latter integrate 
to zero for all cycles on the torus by Stokes' theorem). 
We then employ de Rham's theorem 3 -^, stating that there 
are exactly as many independent holomorphic 1-forms to 
integrate upon the surface as independent 1-cycles to in- 
tegrate along the surface. For the torus the cohomology, 
like the homology, is two-dimensional, i.e. any three (or 
more) holomorphic 1-forms on the torus are linearly de- 
pendent up to an exact form. 



B. Picard-Fuchs equation 

As a result, there must exist a linear combination of 
1-forms {\"(u), A'('u), X(u)} which is an exact form, here 
primes denote derivatives w.r.t. u. This combination 
may be found by allowing for (^-dependent) coefficients 
in front of the three 1-forms and looking for an exact 
form d z [P 2 (z)z~ 1/2 (z 2 + 2uz + 1)~ 1/2 ], where P 2 (z) is a 
second degree polynomial with ^-dependent coefficients. 
Matching coefficients for powers of z leads to 5 equations 
for 6 unknown parameters, determining the sought com- 
bination up to an overall multiplicative factor. This way 
one finds that the operator C = (u 2 — 1)<9 2 + 1/4 acts on 
X(u) as 




FIG. 8: (Color online) Cycles 70 and 71 on the complex 
z-plane for u = —0.9. Notice that cycle 71 crosses twice 
the two cuts from first branch (solid blue line) to second 
branch (dashed red line) and back. 



C\{u) 



d 

dz 



1 

2 z 



1 / 2 (z 2 + 2uz + l) 1 / 2 



(21) 



It follows from Stokes' theorem and the exactness of 
C\(u) that CSj(u) = since jj is a cycle on the torus. 
Thus Sj(u) satisfies the linear second order ODEP^ 



( u i-i)s>;(u) + \s J (u) = o. 



(22) 



This is an example of the Picard-Fuchs equatio n 1 32 * 33 ! (see 
Ref. [34 for a review). Exactly this equation appears 
extensively in the context of Seiberg-Witten theory. 

Inspecting the coefficient in front of the highest deriva- 
tive, one notices that equation (22) has regular singular 
points at u = 00 and u = =pl, where the torus degener- 
ates into a sphere, Fig. [6] Changing variable to u 2 , this 
equation may be brought to the standard hypergeometric 
forrrl 3 -^. In the domain | arg(l — u 2 )\ < tt it admits two 
linearly independent solutions of the form Fo(u 2 ) and 
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uFi(u 2 ), where 

F (u 2 ) 
Fi(u 2 ) 



2F1 
iF\ 



1 1 1 

T"4 ; 2 ; 
1 13 
"4' + 4 ; 2 ; 



(23) 
(24) 



These solutions form a basis out of which Sj(u) (and 
indeed any period of (|14[)) must be composed 



S (u) 
Si(u) 



CooFoiu^ + CoiuF^u 2 ), 
CwFoiu^+CnuF^u 2 ). 



(25) 
(26) 



To find coefficients Cjk, j,k = 0,1 appropriate for the 
action cycles jj one needs to evaluate the periods at 
one specific value of u. Employing the fact that the hy- 
pergeometric functions (23-24) are normalized and an- 
alytic at u = 0, i.e. = 1 + 0(u 2 ), one notices that 
Sj(u) — Cjo + uCji + G(u 2 ). Thus to identify Cjk we ex- 
pand Sj(u) to first order in u and evaluate the integrals 
at u = 0. The corresponding cycles in the z-plane are 
shown in Fig. [9] and explicit calculation yields 



e -i7T/2 Ci 



Coo — e ' 010 
C01 = e + " /2 C u 



87r- 1 / 2 r(3/4) 2 , 
7r-V2 r (l/4) 2 . 



(27) 
(28) 



The relations between Cou and C\u are n °t accidental. 
They originate from the fact that for u = the turning 
points are ±2 and so the cycle 71 transforms into 70 by 
substitution z' = e~ i7T z, Fig. [9} Together with Eqs. (|25|, 



(26) these relations imply global symmetry between the 



two periods 



S (u) = e- i7r/2 S 1 (e i7T u). 



(29) 



Equations (23)-([28j) fully determine the two actions 
So,i(u) through the hypergeometric functions^. One 
should now relate them to physical observables. 



1.5 


Im z 


y^"^ iw 

/V--' 
if u 






-1.5 ; -l.C 1 -0.5 
1 \ 


c.5 1 :.c 


Rez 






**'••...:'€ 

-1.5 





FIG. 9: The two cycles 70,1 for u = 0. Here 71 may be 
mapped to 70 by rotating 180°. 



C. Structure of Sj(u) near u = — 1 

To this end we consider the structure of Sj(u) in the 
neighborhood of u = — 1. As noted earlier, the cycle 70 = 
So contracts to a point as u —> —1 and therefore Sq( — 1) = 
by Cauchy's theorem. By contrast, S±(— 1) remains 
finite. Moreover, while 5o is analytic near u = —1, it 
turns out that 5i is not. To see this, choose some u > — 1 

and allow u to wind around —1 (i.e. (u + 1) — >> (u + 

l)e 2?r *). Since u w — 1 the roots z± in (15) are of the form 



z± = — 1 ± iyj2(u + 1) we see that this transformation 
exchanges these branch points via a counter-clockwise 
half-turn; the branch cut in effect rotates by 180°. For 
the cycle Jo, which encloses the turning cut, this has no 
effect: the cut turns within it. Not so for Si: as the 
cut rotates, one must allow Si to continuously deform if 
Si is never to intersect the branch points. The overall 



effect is shown in Fig. [10] The effect of this monodromy 
transformation is to produce a new cycle S[. Thus, while 
we have returned to the initial value of u, the period 
Si(u) (unlike Sq(u)) does not return to its original value 
and so Si(u) cannot be analytic near u = — 1. 

These facts are consistent, of course, with the origin 
of the integrals as the classical and instanton actions. 
At u —> —1, the classically allowed region collapses and 
p(0) — >• 0, so the classical action at the bottom of the co- 
sine potential approaches that of the harmonic oscillator 
Sq(u) ex (1 + u) (indeed the classical period T ex d u So is 
a constant). For the instanton trajectory 71 the action 
Si does not vanish. Moreover as u — >• — 1 the period on 
the instanton trajectory is logarithmically divergent since 
the trajectory goes to the extrema of the cosine potential, 
Fig.[7| This implies that Si(u) ex const + (l + ^) ln(l + u). 




FIG. 10: (Color online). Monodromy transformation 
(u + 1) — >> (u + l)e 27 ™ rotates the branch cut between 
[2_,z + ] by 180° counter-clockwise. This changes the cycle 
Si — » S'i — Si — So along with it. 

In fact, more can be said. Under monodromy trans- 
formation basis cycle S[ relates to the original basis as 
S'i = Si — So (as may be seen by counting intersec- 
tions of cycles or by moving onto the torus). Thus 
(So, Si) — ^ (Sp, Si — So). From the decomposition of 70, 71 
noted in (20) it follows that the Sj(u) must transform as 



So(u] 
Si(u] 



-2 



'So(uy 



= M_ 



S (u) 
Si(u) 



(30) 
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where we have introduced the monodromy matrix M_i 
of the actions near u = — 1. Since this variation of Si 
occurs for every such monodromy near u — — 1, Si must 
have a component which depends logarithmically on l+u. 
Indeed, In (1 + u) increases by 2iri under the monodromy 
and since Si changes by —2 So it must have the following 
functional form 



SiO) = QiO) + - S (u) ln(l + u) , 



(31) 



where Q\(u) and Sq{u) are analytic functions of (1 + u). 
As an immediate corollary, one can use the relation 



(29) between So and Si to find the structure of the so- 
lution near u = +1. Then the functional form of So{u) 
near u = +1 is So(u) = Qo(u) —iSi(u) ln(l — u)/tv, where 
Qo( u ) = —iQi(-u) and Si(^) = iSo(-u) are analytic 
functions of (1 — u). The corresponding monodromy ma- 
trix is 



Mi 



1 2 
1 



(32) 



While the structure of the periods near u = ±1 has 
been shown through geometric reasoning, it may be also 
found directly by looking for solutions of the Picard- 
Fuchs equation (22 ) as power series in (l±u). Such a pro- 



cedure along with the demand of a constant Wronskian 
leads to a realization that one of the two solutions must 
include (1 ± u) ln(l =b u) terms along with the iterative 
sequence for finding the coefficients of the polynomials. 



This allows for direct verification of Eq. (31 



levels e m = 2au m . To first order one finds for So(u) and 

Qi(u) 

S (u) = 2tt(u + 1), (35) 
Qi (u) = m - - {u + 1) In (32e) , (36) 

7T 

implying e m = —2a + 2a 1 / 2 (m + 1/2). As a result the 
pressure ([5| of a monovalent gas is 

P = -eEoeo = 2k B Tf - ^/k B TeE f. (37) 

The two terms here are the pressure of the ideal gas with 
the fugacity / and the mean-field Debye-Hueckel inter- 
action correction respectively 2 



The instanton action, Eq. (|31|), at quantized u m is 

L/2 



1 \ ( m + 1 
5 1 ( Um ) = 16* + 2i(m+-jln(^r /2 



(38) 



where the linear term in Q\{u) has been absorbed into 
the logarithm. The Gamow formula (34) leads to 



(Ae) m = 2a(Au) m = 2a 



U c ia 1 / 2 S 1 (u m )/2 

i\\fa 



— — ( 326 \ m+1/2 -8a 1 / 2 + (m/2+3/4)lna ( n q) 

"tt \m + l/2j ' [ } 

This coincides wi th the k nown asymptotic results for the 
Mathieu equatiorPESMl. 



D. Semiclassical results 



E. Neighborhood of u — oo 



We now seek semiclassical results for the sequence of 
low-energy bands terminated at u = — 1. We shall inter- 
pret the period Sq(u) which is analytic around u = — 1 
as a classical action. The latter should be quantized ac- 
cording to the Bohr-Sommerfeld rule to determine the 
normalized energies u m of the bands 



S (u m ) = 27ra" 1/2 (ra + 1/2) , m = 0, 1, 



(33) 



(we shall not discuss the origin of the Maslov index 1/2 
here). The second non-analytic period Si(u) is identified 
as the instanton action, which determines the bandwidth 
(Au) m according to Gamow's formula 



Tlx a 



ia 1 / 2 S 1 (u rn )/2 



(34) 



where uj = 2 is the classical frequency for the 
Hamiltonian (12). The monodromy of u around —1, 
Eq. (30), carries over to the bandwidth as a factor of 
e (i/2)a (-250 (n^)) Then the Bohr-Sommerfeld quanti- 



zation (33) is also a condition for the bandwidth to be 



invariant with respect to monodromies. 

To illustrate these results we expand the periods 
Eqs. (25)-(28) near u = — 1 to find the physical energy 



For completeness we also consider the behavior of the 
actions at high energy. In the limit u —> oo the Picard- 
Fuchs equation ([22]) is of the form u 2 S"(u) + S(u)/4 = 
0. Seeking a solution in the form S = u r , one finds 
r(r — 1) + 1/4 = (r — 1/2) 2 = and thus there must be 
two independent solutions with the leading behavior u 1 / 2 
and u 1 / 2 ln(u). So the two periods should be of the form 



Si(u) = u 1/2 [Vi(u) + Wi(u) \nu] 



(40) 



where W$, V{ are analytic functions of 1/u. To find these 
functions one needs to notice that while the continuation 
to infinity for S± is unambiguous, the result obtained 
for So depends on whether the path to infinity passes 
above or below u = 1. This is due to the fact that So 
exhibits nontrivial monodromy around u = 1, Eq. (32). 



In other words, whether u goes to infinity below or above 
the real axis determines which of the two turning points 
z± goes to zero or infinity. Since these are also branching 
points for the torus, the path of analytic continuation 
determines how the cycles on the torus are carried along 
in the process. 

Thus looking for the asymptotic behavior of the peri- 
ods ([25|-([28]) at u -» oo ± iO, one findl^ 
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V (u)=iirW 1 (u)TV 1 (u), (41) 
W (u)=tW 1 (u), (42) 

Vi (u) = Ais/2 [In (e 2 /8) + 2/u] , (43) 

Wi(u) = -4i\/2 [1 - (4zi)" 2 ] , (44) 

to leading corrections in 1/u. Since S'o(w) ± <Si(u) = 
iirWi^u 1 / 2 , from here one may readily show that under 
the monodromy u — > ue 2nl the two actions transform 
with the following monodromy matrices 



M, 



oo— iO 



-3 2 
-2 1 



oo+iO 



(45) 



One may check that the three monodromy matrices sat- 
isfy 

= Mi • M_i, = M i • Mi, (46) 

as expected^ winding around in large counter- 
clockwise circle is the same as winding -1 and 1 sequen- 
tially counterclockwise. 

From Eqs. (40)-(44) one finds the unique non-singular 
period at u — » oo ± iO to be given by So(u) db S±(u) = 
—iiiWiiuju 1 / 2 . As discussed above, it must be identified 
with the classical action and subject to Bohr-Sommerfeld 
quantization (So(u m ) db Si(u m ))/2 = 2ira~ 1 ^ 2 m. This 
leads to u rn « m 2 /2a and thus e m = 2cm m = m 2 , as 
expected for the high energy spectrum. 



V. DIVALENT (2,1) GAS 

The divalent (2,1) gas is the simplest case where the 
Hamiltonian (|8| is non-Hermitian. Employing complex 
variable z = e %e and normalized energy u = 2e/3a, it 
takes the form 



P 



(47) 



Similarly to Eq. (13) this defines a family of complex 
algebraic curves 



f 



T{p, z) = 2p 2 z - (z 3 + 3uz + 2) = 0. (48) 



The map p — p(z) is locally holomorphic away from 
the zeros zo,z± (see Fig. 11). At these three branch- 
ing points as well as at the singularity at z = the 
function p(z) is locally double-valued and behaves as 

respectively. 
" the 



p 



(z — Zj) 1 ! 2 , j = 0, ± and p ~ z 



-1/2 



IV 



In contrast to the monovalent (1,1) case, Sec. 
function p(z) is single- valued at z ~ oo where it goes as 
p ~ z, so no branch cut extends to z = oo. Nevertheless 
there are again four branching points. To construct the 
Riemann sphere we draw two branch cuts: one between 
[0, Zq] and the other between [z+, zJ\. The resulting Rie- 
mann surface is again g = 1 torus, analogous to Fig. [5] 



Its moduli space u contains four singular points u = 
— l,e ±l7r / 3 and u = oo, where the torus degenerates into 
the sphere. (There were only three such points in the 
(1,1) case.) For u = — 1 the branching points z± coalesce, 
while for u = e ±27r / 3 the branching point zq collides with 
z±, correspondingly. As u — > +oo, the branching point 
zq approaches z = 0, while z± —> ±ioc. 




FIG. 11: (Color online) Complex z-plane with two branch 
cuts, shown in gray, (a) Three integration cycles 70,71,72 
are displayed for u = 0. (b) The instanton cycle 
T = —71 + 72. The solid blue (dashed red) lines denote 
parts of the cycles going over the first (second) branch. 

The action integrals are again defined as Sj — §^ A, 

where the 1-form \{u) = p(z)dz/iz is holomorphic on the 
torus. In general the counterparts of the turning points 
in the complex #-plane are not real. This makes it more 
convenient to discuss the action cycles jj in the z-plane. 
With three turning points zo,z±, it is convenient to take 
three paths of integration 70,71,72, depicted in Fig. [TT] 
In terms of the two basic cycles on the torus #o, Fig7j6| 
the three paths are given by 



70 = ^0 : 



7i = -^1 + ^0 , 72 = #1 • 



(49) 



One may notice that 70 — 71 — 72 = 0, and thus 5o = 
Si + S2. This is a consequence of the Riemann- Roch 
theorem^, stating that on a Riemann surface of genus 
1 there are only two independent closed cycles. From 
de Rham's theorem^ one infers that there are exactly 
two independent 1-forms. Therefore the three forms 
{\"(u) : \'(u) : X(u)} are linearly dependent up to an ex- 
act form. Following the root outlined in Sec. IV B| (where 
P2(z) is replaced with Ps(z) - polynomial of degree 3), 
one obtains the Picard-Fuchs equation 



(u 3 + l)S'j(u) + -S j (u) 



0. 



(50) 



In agreement with the above discussion, there are regu- 
lar singular points at the third roots of negative unity, 
i.e. u — — l,e ± * 7r / 3 where the coefficient in front of the 
highest derivative goes to zero, and at u = 00. Two lin- 
early independent solutions Fq(u 3 ) and uFi(u 3 ) of this 
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second-order ODE are given in terms of the hypergeo- 
metric functions 



F\(v :i ) = 2 f, (+^.+^;- w£ 



(51) 
(52) 



In this basis the three periods Sj(u), where j = 0,1,2, 
are given by 



Sj(u) = C j0 F (u 3 ) + CjmF^u 3 ) . 



(53) 



Since the hypergeometric functions Fj(u 3 —> 0) = 1 + 
0(u 3 ), one notices that Sj(u) = C j0 + wQi + 0(^ 3 ), 
as u — >• 0. One can thus find constants Cjk by explicit 
evaluation of the actions at u = 0, i.e. Cjo = Sj(0) and 
Cji = Sj(0). The corresponding integration paths are 
shown in Fig. [TT] and straightforward integration yields: 



Coo = C 10 e™/ 3 = C 20 e-"/ 3 = 
Coi = Cue""/ 3 = C 2 ie" i/3 = 



2 n / 6 37r 3 / 2 

w 

3 1/2 r(|)r(|) 

2 11 /6 7r l/2 ' 



(54) 



(55) 



These relations along with Eq. (53) imply the three- fold 



symmetry between the actions, cf. Eq. (29) 



S (u) = e^Si ( 



e 2i */\ 



(56) 

Now one needs to connect the periods ([53]) with the 
quantum spectrum. We start by discussing the real 
branch of the spectrum terminating at the singular point 
u = — 1, Fig. [I] As u —> — 1, the two branching points 
z± coalesce. As a result 70 cycle degenerates to a point, 
leading to Sq(u — >• —1) — >• 0, while Si^ remain finite 
and actually turn out to be non-analytic. This can be 
seen by considering the monodromy for a windi ng of u 
around -1, i.e. (u + 1) (u + l)e 27ri (cf. Sec. |IVC). 



Such a transformation exchanges branching points z± by 
a counter-clockwise 180° -rotation. This leaves the cycle 
^0 = 7o, which encloses these two points, unchanged. On 
the other hand, the cycle Si picks up a contribution of 
— Sq: S[ = 5i — 5q. Thus 7^2, Eq. (49), pick up a contri- 



bution of ±5q. As a result, for every monodromy cycle, 
5i,2 pick up a contribution of ±So, so locally they are of 
the form 

Si, 2 (u) = Qi, 2 (u) =F f5oW ln(l + u) , (57) 

Z7T 

where Qi, 2(^-0 and So(u) are analytic functions of (l + u) 
(moreover Qi +Q2 = So, cf. Eq. ([49])) . This allows us to 
identify the period 5 (u) = (vW/2)(l + u) + ^) 2 ) 

as the classical action, while the instanton action is a 
combination of the two non-analytic periods Si, 2 (u). 

The corresponding monodromy matrix M_i in e.g. ba- 
sis (So, Si) (since S2 — Sq — Si is linearly dependent) is 



'SoW 



^o(«) 
^l(«) 



= M_i 



So(u' 
Si(w 



(58) 



Employing Eqs. (|49j) , ([56]), one finds that at the singu- 
lar point e i7r/3 (e" 27r/3 ) the period 5i(iz) (S 2 (^)) is non- 
singular and goes to zero. It should be thus identified 
with the classical actions for the branch of the spectrum 
terminating at the respective singular point, Fig. [I] A 
combination of the remaining two actions So and S 2 (Si) 
form the corresponding instanton. The respective mon- 
odromy matrices (again in the basis (So, Si)) are found 
as 



/3 



1 -1 

1 



k/3 



2 -1 
1 



(59) 




FIG. 12: (Color online) Narrow energy bands (red dots) in 
the upper half-plane of complex energy u for a = 200, cf. 
Fig. [3^1. ImSo(u) = along the real axis, where the small 
lines mark ReS (u) = 2na~ 1/2 (m +1/2). The line 
Im Si (u) — emerges from u — e t7T ^ 3 and intersects the 
real axis at u ~ 0.96. The small perpendicular lines mark 
ReSi(u) = 2ita- 1/2 (m+ 1/2). 

To find positions of the bands along the three branches 
of the spectrum, terminating at the three singular points 
u = — l,e ±27r / 3 , one employs Bohr-Sommerfeld quantiza- 
tion for the proper classical action Sj(u) with j = 0, 1, 2, 
correspondingly: 



Sj(u$) = 27ra- 1/2 (m + 1/2), m = 0, 1, 



(60) 



Figure [l2] shows the lines Im So(u) — and Im Si(u) — 
intersected with the set of lines ReSj(u) = 27ra -1 / 2 (m + 
1/2). The numerically computed spectrum sits right at 

the semiclassical complex energies utf) . The excellent 
agreement holds all the way up to the point u ~ 0.96, 
where all three periods Sj happen to be purely real. Be- 
yond this point the semiclassical approximation seems to 
break down, which manifests in e.g. appearance of wide 
Bloch bands. Expanding So(u) near u = — 1, one finds 

for the energy levels e m = 3u$a/2 in the semiclassi- 



cal approximation e m « —3a/2 + y/6a(m + 1/2). The 
corresponding pressure ^ P = —eEoeo consists of the 
two contributions: the ideal (2, 1) gas and the mean-field 
Debye-Hueckel interaction correction. 

Taking into account that there is no physical dif- 
ference between Si and S 2 and that the monodromy 
around u = —1 in Eq. (57) should leave the band- 
width in Gamow's formula (34) invariant (i.e. it adds 
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matrix 



FIG. 13: (Color online) Analytic (numerical) results for the 
logarithm of the bandwidth of the lowest band, m(Ae)o, 
versus square-root of the charge concentration, y/a, with 
(1,1) as dotted (circles), (2, 1) dashed (diamonds) and 
(3, 1) as solid line (stars). 



a factor of exp{(i/2)a 1 / 2 (— 2Sq(u$))}), one identifies 
the instanton cycle with T = —71 + 72, Fig. [TTJ i.e. 
Si nst (u) = —Si(u) + S2(u). This can be also found by 
inspecting the cycles in figure [TT] one sees that the com- 
bined T = —71 + 72 cycle connects z± turning points 
through the "classically forbidden region", similarly to 
71 instanton cycle in (1,1) case, cf. Fig. [8] Expanding 

Si, 2 (u) actions near u = — 1 and substituting u£) from 
the Bohr-Sommerfeld quantization ( |60| ) with j = 0, one 
finds for the Bloch bandwidths of the central spectral 
branch, cf. Eq. (34) with uj = \/6, 



(Ac) 



-a(Au) m 
2V6 ( 36\/6e 



ra+1/2 



(61) 

-3V6a+(m/2+3/4) In a 



m 



1/2, 



Of special interest is the bandwidth of the lowest energy 
band, due to its direct relation to the transport barrier 
of the ion channel, Sec. [Ill Setting m = yields 



(Ae) ~34.14a 3 / 4 e- 7 - 35 ^. 



(62) 



This is in very good agreement with the numerical sim- 



ulations, Fig. 13 



Finally we focus on the behavior at w = 00. The 
Picard-Fuchs equation is of the form u 3 S ff + uS/4 = 0. 
Searching for a solution of the form S(u) = u r leads 
to (r — 1/2) 2 = 0, signifying two independent solutions 
with the leading asymptotic u 1 ! 2 and ?/ 1//2 ln('u). Upon 



the monodromy transformation u 



3 27Ti the first of 



these solutions changes sign, while the second along with 
the sign change picks up a contribution from the first 
one. Considering asymptotics of Si^(u), Eq. (53), at 
u — » +00, one finds the following SX(2, Z) monodromy 



-1 
3 -1 



One can check that 

= M^/s • M_i • M e 



t/3 



(63) 



(64) 



as it should be: winding once around in a large coun- 
terclockwise rotation is identical to winding counterclock- 
wise in sequence around the other three singular points. 



VI. TRIVALENT (3,1) GAS 

The trivalent (3,1) Hamiltonian with the normalized 
energy u is 



u = p 



(65) 



It gives a family of algebraic curves 

£ u : J r (p,z) = 3p 2 z-(z 4 + 4uz + 3) = (66) 

over complex (z,p). They are nonsingular if u A 7^ 1, and 
so T(p^z) implicitly defines a locally holomorphic map 
p = p(z) almost everywhere on (p,z). In this case there 
are six square-root branching points at z = 0, 00 and at 
the four turning points, i.e. four roots of p 2 (z) = 0. 

Hence, while £ u is a doubly-branched cover of the Rie- 
mann sphere, three cuts (instead of two as in the genus- 1 
case) are required per branch. After opening up cuts and 
identifying edges under analytic continuation, this leads 



to a double torus, i.e. a sphere with two handles, Fig.[14fi. 
Unlike the mono- or di-valent cases, the trivalent channel 
gives a family of genus- 2 Riemann surfaces. The excep- 
tional u 4 = 1 cases make £ u singular at (p, z) — (0, —u), 
due to collision of two turning points, Fig. [T4) 3. So the 
double torus degenerates into a simple torus with two 
points identified (a singular surface of genus 1). 




FIG. 14: (Color online) (a) Double torus curve &u 

for u 4 / 1, 

having four basic cycles, (b) When u A = 1 the g = 2 torus 
degenerates into a singular g — 1 surface. This makes one 
of the basic cycles to pass through the singularity, and 
renders another cycle contractible to a point. 



As in the genus- 1 cases, the actions can be understood 
as integrals Sj = §^ A of the holomorphic action 1-form 
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X(u) = p(z)(dz/iz) upon these Riemann surfaces. Owing 
to the four turning points, there will be four such cycles 
7j with j = 0, 1, 2, 3. These are chosen as in the divalent 
case, with the inner arcs of each being taken to start 
on the principal branch. They are shown for u = in 
Fig. [T5fc . The ^-dependence of these periods is governed 
by the Picard-Fuchs equation. 



lutions of the form u k Fk(u 4 ), where k = 0, 1, 2, 3 and 





Im z 














72 f / 




IjRez 








T7o 
















(a) 


(b) 



FIG. 15: (Color online) The Riemann surface is doubly 
branched with a total of three cuts, shown in gray. The 
four cycles with j = 0, 1, 2, 3, along with the instanton 
cycle T (defined for later reference) are displayed for 
u = 0. The solid blue (dashed red) lines denote parts of 
the cycles going over the first (second) branch. 



As the double torus is genus-2, there are four in- 
dependent cycles (as opposed to two for genus-1). 
So the homology — and so too, as argued before, the 
cohomology — is not two- but /b^r- dimensional: any five 
holomorphic 1-forms on the double torus are linearly de- 
pendent up to an exact form. Thus X(u) and its first four 
derivatives can be used to produce an exact form; this is 
done by finding coefficients in a polynomial entering the 
exact form, as discussed in Sec. |IVB| Stokes' theorem 
implies that S(u) = §^ \{u) must satisfy a 4th-order lin- 
ear ODE in u, i.e. Picard-Fuchs equation which in the 
present case takes the form 



(u 4 - 1)S^ + 8u 3 S^ + 



217 

"18" 



i 2 S" 



uS' 



65 
144 



S 



= 0. 
(67) 



It has regular singular points at fourth roots of 1, i.e. 
u G {±l,±i} and at u = oo. By changing variable to 
u 4 , one can cast the Picard-Fuchs equation as a gen- 
eralized hypergeometric equation. In the cut domain 
| arg(l — u A )\ < 7r it has four linearly independent so- 



F (u 4 ) 



4F3 



F 2 (u 4 ) 



F 3 (u 4 ) 



4^3 



4^3 



1 5 13 1 


1 3 


8 ' 24 ' 24 ' 4 ' 


2' 4 


1 11 in 1 

1 11 iy 1 





8' 24' 24 ' 2' 


4' 4 


3 17 25 3 


5 3 


8' 24' 24 ' 4' 


4'2 


5 23 31 5 


3 7 


8' 24' 24 ' 4' 


2' 4 



(68) 



(69) 



(70) 



(71) 



are generalized hypergeometric series. Note that the pa- 
rameters of each 4,F^{{ai}] u A ) satisfy — Xl a * = 
1; such hypergeometric series are known as one-balanced 
or SaalschutziarP^. 

Writing the actions in this basis as 



Sj(u) 



3 

E 

fc=0 



C jk u k F k (u 4 



(72) 



we note that Sj(u) = Yll=o CjkU k -\-0(u 4 ) (as generalized 
hypergeometric functions are unity at zero and analytic 
nearby). We expand each Sj(u) up to u 3 around u = 

to obtain 



(73) 
(74) 
(75) 
(76) 



and evaluate the resulting integrals, Fig. [L5[ 
the {Cjfc}^l For e.g. So this brings 

Coo = +2 7/2 • 3- 9/8 7r" 1 / 2 r(5/8)r(7/8) 
C01 = -2" 1/2 • 3- 7/8 7r- 1/2 r(l/8)r(3/8) 
C 02 = -2" 5/2 • 3- 13/8 7r- 1/2 r(l/8)r(3/8) 
C 03 = +7 • 2" 1/2 • 3- 27/8 7r- 1/2 r(5/8)r(7/8) 



When u = the turning points satisfy z 4 + 3 = and so 
they lie on a certain circle in the complex plane. Hence 
7j and 7^+1 are only different by tt/2 rotation, Fig. 15 1. 
As a result, we find the four- fold symmetry relations 

S (u) = e?Si(e-?ti) = e^S 2 (e- 7Ti u) = e-^S^e^u) 

(77) 

for u in the cut domain | arg(l — u 4 )\ < tt. 

We now consider the periods in the neighborhood of 
u = — 1. As before, the cycle 70 becomes contractible 
to a point as u — » —1 and therefore Sq{— 1) = by 
Cauchy's theorem. The other three actions remain fi- 
nite, but Si and 63 are non- analytic. This can be seen 
by considering the monodromy around u = — 1. As in 
the genus-1 cases, the shrinking branch cut near z = 1 
makes a half- turn. Examining the action cycles, it is 
only 71 and 73 that intersect the cut rotating under the 
monodromy within the 70 cycle. Hence it is these two 
cycles that change under monodromy and thus have log- 
arithmic non-analyticity near u = — 1. More precisely, 
(Si, S3) — >> (Si + So, S3 — So) under the monodromy and 
so these actions are of the form 



Si, 3 (^) = Qi, 3 W T 7T S o(u) m(l + u) , 

Z7T 



(78) 
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where Q 1,3(11) as well as So(u) and 82(11) are analytic 
near u = — 1. Since + 63 is seen to be invariant 

under the monodromy, there are a total of three inde- 
pendent periods which have trivial monodromy around 
u = — 1. This is again supported by considering series 
solutions of the Picard-Fuchs equation (67) near u = — 1. 



This way one finds three regular solutions with leading 
behavior (1 + (1 -f iz) 1 , (1 + u) 2 along with an irreg- 
ular solution with the leading behavior (l-\-u)\n(l-\-u). 
For reasons of space we omit the corresponding 4x4 
monodromy matrix. 

Althou gh analytical facts about 4F3 series are sparse 
(se d 35 | 3S | for relevant discussion), there are simple consis- 
tency checks which our solutions (72) must pass. First 
the vanishing of the classical action Sq(u) at w = 
implies the identity 



1 



^C ok (-l) k F k (l) = 



(79) 



k=0 



for the hypergeometric functions given above. In addi- 
tion, from inspection of the Hamiltonian (65 ), one notices 
that the classical frequency near u = — 1 is oj = VS. This 
implies S' ( — l) = ^(2tt/lj) and thus 



k=0 



>7T 



(80) 



Being checked numerically, both hold up to 10 -16 . 

Now we turn to the analysis of the spectrum of the 
Hamiltonian (65) at large a. There are three spectral 
branches terminating at the singular points u = — l,±z, 
Fig. [3)3 (notice that the fourth point u = 1 lies in the mid- 
dle of the spectrum and does not have an obvious semi- 
classical interpretation). To determine positions of the 
bands we quantize the corresponding actions j — 0, 1,3 
(but not j = 2, the latter is responsible for the period 
vanishing at u = 1) according to the Bohr-Sommerfeld 
rule: 



Sj(u$) = 2^a- i r\m + 1/2), m = 0, 1, j = 0, 



-1/2 



Figure 
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shows the semiclassical energies along with 
numerically found energy bands. One notices the perfect 
agreement between these two for Reu < 1.09. At the 
point u ~ 1.09 all three actions So, 1,3 are purely real 
and the corresponding instanton action (see below) goes 
through zero. Beyond this point energy bands are not ex- 
ponentially narrow and semiclassical approximation may 
not be applicable. Notice that this point is unmistakably 
different from the singular point u = 1. Focusing on the 
real energies at the bottom of the spectrum and expand- 
ing near u = — 1, one finds with the help of identities ( 79 ), 
@ S (ii) = (v / 87r/3)(l + u) + 0(1 + u). The Bohr- 



1,3. 
(81) 



Sommerfeld rule (81) leads to e m = A.u£)a/% = — 4a/3 + 
2 v / 2a 1/2 (m + 1/2). Employing Eq. ([H), this yields the 



The two terms here are the ideal gas pressure and the 
mean-field Debye-Hueckel correction respectively. 




FIG. 16: (Color online) Narrow energy bands in the upper 
half-plane of complex energy u for a = 200, cf. Fig. |3]d. 
ImSo(u) = along the real axis, where the small lines 
mark Re S (u) = 2na~ 1/2 (m+ 1/2). The line 
Im Si (u) = emerges from u — i and intersects the real 
axis at u w 1.09. The small perpendicular lines mark 
ReSi(u) = 27va~ 1 / 2 (m + 1/2); red dots, numerically 
computed narrow bands. 



Let us now focus on the width of the Bloch bands near 
u = — 1. This requires to identify a cycle corresponding 
to the instanton action. Guided by the cosine potential 
example, cf. Fig. [7| we take the corresponding cycle as 
connecting the turning points of the classical action So 
through the "classically forbidden region" . This suggests 
cycle T shown in Fig. [T5] o ? which is essentially of the 
same form as 71 instanton cycle in (1,1) case. One can 
see that T = 73 — 72 — 71 by considering intersections 
of these cycles. Upon the monodromy transformation 
around u = — 1 the instanton action thus acquires a con- 
tribution — 2So(u), Eq. (78), which leaves the bandwidth 



invariant thanks to Bohr-Sommerfeld quantization (81). 
The resulting instanton action is 



Sinst(u) = Qinst(u) + ~S (u) ln(l + u) , (82) 

7T 

where Qi ns t = Q3 — S2 — Qi is the regular part of 



Sinst(u), cf. Eq. (78). To first order in (1 + u) this is 
Qinst(um) ~ 14.12z — 6.7H • (1 + u), where e.g. the lead- 
ing term originates from 

Qinst(~ 1) = $>inst(~ 1) 
3=0 

Then, for u^ni along the real ?i-axis satisfying Bohr- 
Sommerfeld quantization, Gamow's formula yields for 
the bandwidth 



(A e ,_f,A„,_f.J^,-™ 
4V2 ( 6.35 \ m+1/2 



pressure of the trivalent Coulomb gas as P = %ol — y2a. 



1/2 



(83) 

7.06Va+(m/2+3/4) In a 
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The width of the lowest band (Ae)o is compared with the 
numerical results in Figfl3l As in the earlier cases the 
two results are in strong accorcP^. 

For completeness we address the u = oo behavior. 
For large u the Picard-Fuchs equation is of the form 
u 4 S^ +8u 3 S^ + 217u 2 S" / 18 + uS' + 655/144 = 0. The 
trial S(u) = u r brings four independent solutions with 
leading asymptotic {u 1 / 2 , u 1 ! 2 ln(V), u~ 5 / 6 , u~ 13 / 6 }. The 
former two are familiar from the genus- 1 cases, but the 
latter two are novel to the genus-2 case. The fractional 
powers ex 1/6 may seem unexpected, given the four- fold 
symmetries of the periods. However, this symmetry is 
manifest at the level of cycles at u — 0, where four turn- 
ing points are equally spaced on a circle in the complex 
z-plane. By contrast, as ii — ^ oo, the turning points 
must satisfy either z 3 ~ — u or 1/z ~ — u, thus only 
three of the four turning points tend towards infinity and 
one towards zero. This leads to the three- fold exchange 
of actions upon monodromy around u = oo. Thus the 
u r behavior of the periods with r = — integer/(2 * 3) is 
exactly what is needed to construct a proper SX(4, Z) 
monodromy matrix. 



VII. HIGHER VALENCE GASES 

Here we briefly summarize our current state of under- 
standing of the higher valence (4, 1) and (3, 2) gases. The 
corresponding Hamiltonians are 



(4.1) : 

(3.2) : 



u = p 



■ u = p 



z 4 1 
3 + 2z 2 ) ' 



(84) 
(85) 



In both cases there are five turning points in the z-plane 
given by the equation p 2 (z) = 0. The behavior at z = 
and z = oo is somewhat different: for (4, 1) there is a 
branching point at z = 0, but not at z = oo (cf. (2,1) 
problem); while for (3, 2) the opposite is true: there is no 
branching point at z = 0, but there is one at z = oo. In 
either case there are six branching points, which dictate 
three branch cuts. The resulting Riemann surface is the 
double torus, as in (3, 1) case, Fig. 14 In these cases it 
is not degenerate as long as u 5 ^ —1; otherwise two of 
the five turning points collide, leading to a contraction of 
one of the cycles. Therefore one expects five branches of 
the spectrum terminating at u = (— l) 1 / 5 , in agreement 
with Figs. [3f;,d. 

Since the Riemann surfaces are genus-2, there is a lin- 
ear combination of the 1-form X(u) = p(z)dz/iz and its 
four ^-derivatives which sum up to an exact form. There- 
fore any period S = § A must satisfy a 4th-order ODE in 
u. This is found by matching coeffi cients in a polyno- 
mial entering the exact form (see Sec. IV B), yielding the 



Picard-Fuchs equations 

Cb/ 5 — 1 

(4,1): (u 5 + l)S^(u) + — -S^(u) 



(86) 



+ ^ u 3 S" (u) + \ u 2 S'(u) + § uS(u) = 0, 
lb 4 o4 



— 1 

(3,2): (u 5 + l)S^(u) + — -S^(u) 



(87) 



+ ±f u 3 S"(u) + \ u 2 S'{u) + ^ uS(u) = 0. 



While the coefficients seem arbitrary, some features are 
notable. First, changing variable to the equations can 
be brought to the generalized hypergeometric form; one 
finds four independent solutions of the form u k Fk(u 5 ), 
where k = 0,1,2,4 and being a certain 4 F 3 hyper- 
geometric seriePM Notice the absence of a k = 3 solu- 
tion. This can be verified directly from the Picard-Fuchs 
equations, whose leading behavior near u = is given by 
SW(u) -u^S^iu) = 0. Substituting S oc u k , one finds 
fc(fc — l)(fc — 2)(fe — 4) = 0. 

Second, let us focus on the vicinities of fifth roots of 
-1, e.g. on u = -1. Notably both Eqs. @,([87} have 
the same leading behavior b(u + l)S^(u) + 10S&>(u) = 
0, with all other terms are subleading. Looking for a 
solution in the form S(u) ~ (1 + u) s , one finds for the 
8-exponent 5s(s — l) 2 (s — 2) = 0. Therefore in both 
cases there are three analytic solutions with the leading 
behavior (1 + ^)°,(1 + ^) 1 ,(1 + ^) 2 , while the double root 
at s — 1 signifies that the fourth independent solution is 
of the form (1 + u) ln(l + u^. 

This observation indicates non-trivial monodromy ma- 
trix M_i, allowing one to identify the polynomial in front 
of the ln(l + u) with the classical action Sq(u). Being 
quantized according to Bohr-Sommerfeld, the latter de- 
termines the spectrum along the branch terminating at 
u = -1, Figs.^d. 

Finally, we consider the behavior at u — >• oo. By taking 
trial solutions in the form S(u) ~ u r , one obtains 4-th 
order algebraic equations for the exponent r. The four 
roots of these equations are {^,^,— ^} f° r (4, 1) 
case and {^,^,— |,— ^} for (3,2) case. Remarkably, 
there is a double degenerate root at r = 1/2 in both cases, 
leading to the two solutions with the leading asymptotic 



behavior v}l 2 and ^ 1 / 2 ln(^). This was also the case in 
all the examples, considered above. The first of these 
solutions, being quantized, leads to e m = m 2 , expected 
at large energies. The two other roots bring two addi- 
tional solutions with the leading behavior ?x -3//4 , u~ 13 ^ A 
or w -7 / 6 , u~ 17 / 6 for (4,1) and (3,2) cases, correspond- 
ingly. The denominators of these fractional powers may 
be related with the fact that four and three turning points 
go to infinity as u — >• oo in the two respective cases. The 
monodromy transformation interchanges the corre- 
sponding periods (possibly with a sign change). This is 
achieved by having —integer/4 and —integer/ (2 * 3) pow- 
ers in the corresponding solutions. 
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VIII. CONNECTIONS TO SEIBERG-WITTEN 
SOLUTION 

Here we briefly rev iew t he main features of Seiberg- 
Witten (SW) solutioiPEHl which were adopted in our 
calculations^. The original SW construction gives the 
spectrum of a four-dimensional supersymmetric SU(2) 
Yang Mills theory (SYM). Spectrum of the infrared the- 
ory appears to be given by the set of electrically and 
magnetically charged particles (BPS dyons), which are 
different from the fundamental particles of the initial UV 
theory. The latter consists of a vectormultiplet trans- 
forming in the adjoint representation of SU(2), whose 
components are: one complex scalar field 0, pair of Weyl 
fermions (gluini) and a SU(2) gauge field (gluon). In a 
classical UV vacuum <\> aligns along the Cartan generator 
of 5u(2) as ((/)) = (10-3/2, where the complex expectation 
value a parameterizes the manifold of classical vacua. In 
the quantum theory a more convenient coordinate is 



(tr0 2 ) 



(such that in the classical limit u —> 00 one has u ~ a 2 ), 
defining the moduli space of quantum vacua of the theory 
M„. 

Given the expectation value a, one defines the gen- 
erating function (prepotential) J 7 (a) as a logarithm of 
the partition function of the theory, restricted by ((f)) = 
(MJ3/2. It allows to introduce a canonically conjugated 
complex variable 



a D 



dJ r (a) 
da 



(89) 



where one may regard (a, op) as the coordinate and mo- 
mentum on A4 U . The underlying supersymmetry allows 
to argue that a(u) and ajj(u) are holomorphic functions 
on the moduli space, safe possibly for few isolated singu- 
lar points. In the UV limit u — »• 00, one finds a one-loop 
correction of the form 



a D 



la 

7T 



1 + In 



A 2 



(90) 



where A is a dynamical scale. Recall that a ~ ^fu in 
this region. Therefore, when the argument of u changes 
by 27ri, a changes its sign and transforms as a/) ^ 
— od + 2a. This rule can be parameterized using the 
following monodromy matrix in the (a.o,a) basis 



Moo = 



-1 2 
-1 



(91) 



To find the spectrum of the IR theory means to com- 
pute masses of particles which are protected by super- 
symmetry (so called BPS dyons). BPS mass formula 
reads 



where (n e ,n m ) are electric and magnetic charges of a 
dyon respectively, e.g. a monopole has (n e ,n m ) = 
(0, ±1). The above relationship can be understood semi- 
classically (at large u) by evaluating the energy func- 
tional for the UV theory on the electrically and magneti- 
cally charged configurations. The AT = 2 supersymmetry 
guarantees that the very same formula works at strong 
coupling as well. There are special loci in the u plane 
where the masses (92) vanish. One can identify these 



points as singularities for a and an- 

Let us look at the point uq, where the monopole be- 
comes massless cld(uo) = 0. By a conformal transforma- 
tion one may always scale uq = 1. In a vicinity of this 
point an behaves as cld oc (u — 1), thus near this point 
d]j(u) is holomorphic, while a(u) is expected to be sin- 
gular. Performing a one- loop calculation similar to the 
one near u = 00, in the framework of dual theory, one 
obtains a relation similar to ( 90 ) 



ia D a D 
a ~ m — — 

7T A 



(93) 



Recalling that an ~ (u — 1), one finds for the monodromy 
matrix near u = 1, again in (cld, a) basis: 



Mi 



1 

-2 1 



(94) 



From the symmetry considerations one may argue that 
there should be at least one more singularity in addition 
to u = 00 and u = 1. It follows from the fact that if 
a singularity exists at some value of uq there ought to 
be another one at — uq. The Z2 symmetry, which flips 
the sign of u, is a result of breaking of the global U(l) 
symmetry (so-called R-symmetry) of IR action. The lat- 
ter is a remnant of the analogous symmetry in the UV 
theory which is common for gauge theories with an ex- 
tended supersymmetry. It exists on the classical level, 
but is broken by quantum corrections (both perturbative 
and instanton) down to the Z2 for u = (tr</> 2 ). There- 
fore, there are at least three singularities in A4 U , e.g. at 
u = 00 and u = ±1. The third singular point u = — 1 cor- 
responds to a massless dyon of unit electric and magnetic 
charges a(— 1) + od(— 1) = 0. The monodromy matrix 
around it can be computed employing completeness re- 
lation M\M-\ = Mqo in the complex i^-plane. 

The non-trivial realization of the SW construction is 
that complex variables (ao(u), a(u)), with the analytic 
properties deduced above, may be viewed as periods of 
algebraic curves (tori) £ U: defined over the moduli space 
M. u , with respect to some meromorphic differential Xsw- 
The simplest way to parameterize such a curve is 



T(y, x) = y 2 - (x - u)(x - l)(x + 1) = , (95) 



M nejnm (u) = \n e a(u) + n m a D (u)\ , 



(92) 



where x, y are complex. The above equation describes a 
double cover of the x-plane branched over the four points 
x = ±l,u and x = 00. Moreover the cover is singular 
any time two of these points coalesce, i.e at u = ±1, 00, 
as required. Basis in the first (co) homology of £ u (two 
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dimensional in this case) is given by integrals of a one- 
form over one-cycles. We pick the homology basis So, Si, 
Fig. [6l and one- form Asw(^) (SW differential) such that 



a D (u) 



A S w , a(u) 



So 



Si 



Asw • (96) 



To pick a proper SW differential Asw(^)? we recall that 
there are only two linearly independent holomorphic 1- 
forms on the torus up to an exact form. These two forms 
may be chosen as Ai = dx/y and A2 = xdx/y, so Asw = 
Pi( u )y +/M 1 0^5 where f3i j2 (u) are functions of u only. 
The requirement that the period integrals ( [96] ) reproduce 
correct asymptotic behavior of a(u) and ahyu) at u = 1 
and u = 00 (90), (93) allows to determine fii^iu). Finally 
one obtains 



Asw = 



V2 

2k vx 



1 



: dx . 



(97) 



From here one can evaluate the periods ( 96 ) in terms of 



elliptic integrals. They in turn yield the entire informa- 
tion about BPS mass spectrum (92) and the prepotential 
p). 



Close parallels to our calculations are apparent. In fact 
the SW construction, outlined above, essentially mirrors 
the (1, 1) gas calculations. The elliptic curve ( [95]) is iso- 
genic to the torus ( 14 ) and the two SW periods (|96[) are 
directly related to the two action integrals as Sp y dp 
and S\ ~ a + cld- In fact, they may be showrP^^ to 
satisfy exactly the same Picard-Fuchs equation ([22]) as 



our actions. Therefore the two basis solutions (|25|), (|26|), 
expressible through the complete elliptic integrals of the 
first and second kincPSI, are also a basis for SW periods 
ar>(u), a(u). 

An interesting open question is whether our multiva- 
lent examples have analogs in SYM theories. For exam- 
ple, (2, 1) case corresponding to a torus with the residual 
Z3 symmetry in the w-plane, may be related to SU(2) 
theory with several fundamental hypermultiplets added. 
Other examples, leading to g = 2 surfaces with Z4 and 
Z 5 symmetries may be related to certain SU(3) SYM 
theories with matter. 

Another captivating observation is related to the pe- 
culiar structure of the spectra near u « 0.96 in (2, 1) 
gas, u ~ 1.09 in (3, 1) gas, etc. These points are marked 
by the condition Im Si{u)/ So(u) = 0, which is reminis- 
cent of wall crossing phenomena in J\f = 2 theories^. 
It is observed that moduli space M u has domains sepa- 
rated by walls such that when one "crosses" a wall the 
spectrum of the IR theory changes dramatically. For in- 
stance, for the SU(2) theory at small \u\ there are only 
two states in the spectrum: monopole (0, ±1) and dyon 
(±1,=f1). However, at large \u\ these particles can form 
bound states with higher electric charge (n, ±1) for any 
integer n. The wall is given by ImaD(u)/a(u) = 0. 



IX. DISCUSSION OF THE RESULTS 

In this paper we developed semiclassical treatment for 
a family of non-Hermitian VT- symmetric Hamiltonians. 
These Hamiltonians appear upon transfer-matrix map- 
ping of ID classical statistical mechanics of multi-valent 
Coulomb gases onto quantum mechanics. The low-energy 
spectra of the Hamiltonians directly translate into ther- 
modynamic and adiabatic transport coefficients of the 
corresponding Coulomb gases. 

We use methods of algebraic topology, traditionally 
employed in the context of the Seiberg-Witten theory. 
The main advantage of this strategy is that it allows us 
to avoid solving equations of motion and finding classi- 
cal trajectories explicitly. The latter task is rather non- 
trivial (if at all attainable) in the 4D phase space. In- 
stead, we argue that any surface of constant energy is a 
2D Riemann surface with genus g > 1. The action along 
any closed trajectory (not necessarily satisfying equations 
of motion) may be written as an integer valued linear 
combination of 2g basic periods of the surface. The lat- 
ter may be found as solutions of Picard-Fuchs ODE in 
the space of parameters. Finally, relations between ba- 
sic periods and the quantum spectra are established by 
considering special points in the parameter space, where 
the surface degenerates into genus g — I singular surface. 
Consideration of monodromy transformations in a vicin- 
ity of these points allows us to identify classical actions, 
quantized according to Bohr-Sommerfeld, as well as the 
instanton action, which determines the bandwidth. 

Results obtained this way are in excellent agreement 
with numerical simulations in a broad range of param- 
eters. One of the reasons for this success is that the 
method provides with preexponetial factors on the same 
footing with the exponent itself. Another appealing fea- 
ture of the approach is that none of our semiclassical cal- 
culations required the concept of imaginary time. In fact 
"time" (i.e. ID coordinate of the Coulomb gas) does not 
appear at all. In a sense it is substituted by evolution 
in the space of parameters of the Hamiltonian (moduli 
space). We expect the method to be useful in a broad 
class of problems which require instanton calculations in 
complex spaces. 
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